Hydrodynamic theory of de-wetting 
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A prototypical problem in the study of wetting phenomena is that of a solid plunging into or being 
withdrawn from a liquid bath. In the latter, de-wetting case, a critical speed exists above which 
a stationary contact line is no longer sustainable and a liquid film is being deposited on the solid. 
Demonstrating this behavior to be a hydrodynamic instability close to the contact line, we provide 
the first theoretical explanation of a classical prediction due to Derjaguin and Levi: instability 
occurs when the outer, static meniscus approaches the shape corresponding to a perfectly wetting 
fluid. 



The forced wetting or de- wetting of a solid is an impor- 
tant feature of many environmental and industrial flows. 
In typical applications such as painting, coating, or oil 
recovery it is crucial to know whether the solid will be 
covered by a macroscopic film or not. If a fiber is plung- 
ing into a liquid bath to be coated (wetting case), the 
speed can be quite high (m/sec) while maintaining 
a stationary contact line. In the opposite case of with- 
drawal (de- wetting) , a stationary contact line is observed 
only for very low speeds, and a macroscopic film is de- 
posited 0,13 typically at a speed of only a few cm/sec. 
Yet no theoretical explanation for this instability exists, 
or of the fundamental difference between the two cases. 



It is well known [4], |5j that viscous forces become very 
large near a moving contact line, and are controlled only 
by some microscopic cut-off A, for example a slip length 
4, 6]. As a result of the interplay between viscous and 
surface tension forces, the interface is highly curved, and 
the contact line speed U is properly measured by the 
capillary number Ca — Urj/j, where rj is the viscosity of 
the fluid and 7 the surface tension between fluid and gas. 
Owing to this bending the interface angle measured at, 
say, away from the contact line differs [?l IsL I9I Ho| 

significantly from the microscopic angle directly at the 
contact line. 

It was first proposed by Derjaguin and Levi an d 
later reiterated by others 0], that instability occurs if 
this dynamic interface angle reaches zero. Apart from 
the fact that there is no justification for this condition, 
it does not lead to a unique criterion since the angle de- 
pends on the position where it is evaluated. However, it 
was noted experimentally 2] that the interface profile at 
the critical speed corresponds to a static meniscus with 
zero equilibrium contact angle, except in the immediate 
neighborhood of the contact line. Another important set 
of experimental observations [l^j ] was performed with flu- 
ids of different viscosities and equilibrium contact angles 
9 e in a capillary tube: it was found that instability oc- 
curred if Caj9\ exceeds a critical value, strongly pointing 
to a hydrodynamic mechanism for the instability. 

To explain the physical mechanism behind the insta- 
bility, we note that the interface shape h{x) near the 



contact line must have the form 



h(x) = 3Aff(0, £. = x0 e /(3\), 



(1) 



since the cut-off A is the only available length scale. The 
dependence on 9 e was put in for later convenience. As 
expected, the curvature h"(0) = 0g/(3A)if"(O) becomes 
very large at the contact line, since A is in the order of 
nanometers. As the curvature of the local solution has to 
match to a value of order unity in the static outer part of 
the profile, the usual boundary condition for H is one of 
vanishing curvature for large £. In the wetting case this 
leads to an asymptotic solution first found by Voinov 0] 
and an expression for the interface angle usually referred 
to as "Tanner's law" Q. 

However, it is a well-known but little appreciated fact 
[lfij that Voinov's solution fails away from the contact 
line in the de-wetting situation. Instead, the local so- 
lution always retains a positive curvature, and will fail 
to match to an outer solution in the limit of small A. 
At the critical speed, the necessary compromise between 
the inner solution near the contact line and the outer 
static solution has been pushed to the limit: from all 
possible inner solutions, the one with the smallest pos- 
sible curvature is selected. The outer solution, on the 
other hand, realizes the solution with the largest curva- 
ture, which happens to be the one corresponding to zero 
contact angle, in agreement with the criterion of Der- 
jaguin and Levi. Note that this mechanism differs from 
the one proposed by de Gennes [lfij . which is based en- 
tirely on properties of the local solution near the contact 
line. In a companion paper [l7| we show that de Gennes' 
local solution is inconsistent with a class of contact line 
models, such as the slip model considered here. 

To lay down the theory behind the above argument, 
we choose a plate being withdrawn from a liquid bath 
(see FigQJ as a test case, a geometry which was stud- 
ied in detail in ^^|. There it was found numerically 
that there exists a critical capillary number Ca cr , above 
which no static solution exists, in agreement with the ex- 
perimental findings described above. If both 9 and the 
equilibrium contact angle 9 e are small, one can treat the 
problem in the framework of lubrication theory, assum- 
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moving plate 



FIG. 1: A schematic of the setup. At the contact line, h(0) 
= 0, and the slope of the interface is 6 e . The plate is being 
withdrawn at an angle 6. 



ing a parabolic flow profile. This greatly simplifies the 
problem, but without altering its essential structure. 

To relieve the corner singularity at the moving contact 
line, we allow the fluid to slide across the solid surface, 
following the Navier slip law 



du 

u| y =o -U = X Q^\y=o 



(2) 



at the plate, where h(x) is the thickness of the fluid layer. 
The resulting lubrication equation is |l8j 



3Ca 



h 2 + 3Xh 



ti" - ti 



(3) 



where we consistently used the small-angle approxima- 
tion tan(0) ~ 9. All lengths are scaled by the capillary 
length l c = y/j/{pg). 

The left hand side of J3J corresponds to viscous forces, 
diverging as the contact line position h(0) = is ap- 
proached, but weakened by the presence of slip. Close 
to the contact line, viscous forces are balanced by sur- 
face tension forces (first term on the right), resulting in a 
highly curved interface near the contact line. The other 
two terms stem from gravity and only come into play 
at greater distances. We also assume that the angle at 
the contact line h'(0) — 9 e is constant, which then has 
to be the equilibrium contact angle, in order to give the 
right result at vanishing speed. Note that we have used 
the simplest possible contact line model, and no special 
significance should be attached to our choice. Since the 
instability mechanism is hydrodynamic, any local contact 
line model can be incorporated into our description. Far 
away from the interface the surface coincides with the liq- 
uid bath, so the third boundary condition is h'(oo) = 9. 

To disentangle the local contact line behavior from the 
far-field meniscus, we transform equation J2J into simi- 
larity variables l(T]l: 



H 2 +H 



H'" + n 2 {9/9 e - H') 



(4) 



with boundary conditions H(0) — 0,H'(0) — 1, and 
H'(oo) = 9/9 e . The three parameters remaining in the 
problem are the rescaled capillary number 8 = 3Ca/9^, 
the rescaled cut-off /i = X/9 e , and the relative inclination 
angle 9/9 e . 

Since we are interested in the limit of (i being very 
small, we will first look at the equation for fi — 0, as ap- 
propriate close to the contact line. As mentioned above, 
in the limit of large £ this equation behaves very differ- 
ently, depending on whether 8 > or 8 < 0. This can 
be understood from the simplified equation H'" = 8/ H 2 , 
which can be integrated completely [lil . Namely, if 
8 < (wetting), there exists an asymptotic solution for 
£ — > oo j 1 11] - whose curvature vanishes at infinity. 

If on the other hand 8 > (de- wetting) , all solutions 
have a finite curvature at infinity. Solutions which o bey 
our boundary condition H(0) = can be written |l5| 
parametrically as 



f. _ 2 1 ^7rAi(s) 

? ~ P(uAi(s)+13Bi(s)) 

IT ' 

n ~ (uAi(s)+[iBi{s)y 2 



s 6 [sioo[, 



(5) 



where Ai and Bi are Airy functions. The limit £ — » 
corresponds to s — > oo. For large £, the curvature 
becomes fl5| 



^irAi(si) 
where s x solves the equation 



>0, 



aAiisx) + PBi(si) = 0. 



(6) 



(7) 



Since the solution extends to s = oo, s± has to be the 
largest root of (JJJ . 

The constant f3 can be determined by matching (J5J to 
the cutoff region near the contact line, where one finds 
« [3<51n(2 2 / 3 /3 2 /«)] 1/3 in the limit of small £ 
|l5|| . Comparing this to the first order expansion of the 
full equation (JIJ near the contact line 



r expar 

HQ, 



£T'(0 = !-*(!+ MO). 



(8) 



we find (3 2 = 7rexp(l/(3<5))/2 2 / 3 + 0(8). The matching 
described here was investigated in greater detail for the 
case 8 < in 01 ■ We found that higher order corrections 
in 8 were surprisingly weak, and depended only very little 
on the type of cut-off used at the contact line. Thus we 
are confident that the same holds true in the present case. 

At the critical capillary number Ca cr the only remain- 
ing parameter, which is a, can be determined from the 
condition that should be minimal, in other words 
Ai(si) must be maximal among solutions of 0. By 
choosing a — a cr = — (3Bi(s max ) / Ai(s max ) we can in fact 
ensure that Ai assumes its global maximum 0.53566 . . ., 
which occurs for s — s max — —1.0188 . . .. In summary, 
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we find then for the critical solution that minimizes the 
curvature 



(5 1 / 3 exp[-l/(3(5)] 
2V3 7r ( i 4i( w ))2 = 



(9) 



which is now given exclusively in terms of the rescaled 
capillary number 8. 

The solution © has to be matched to an outer solu- 
tion of (0 , valid away from the contact line. In the spirit 
of the classical technique |(| employed for the spreading 
drops, this is a static solution of (J3J), but with an appar- 
ent contact angle different from H'(0) = 1. The general 
solution of J3J with 5 = gives 



H'(£) = 6/6 e - 0o e 



(10) 



so the curvature at the plate is H"(0) = iiQq. The largest 
curvature for which (|10|l still makes sense is fJ,0/0 e , at 
which point the apparent contact angle has reached zero. 

The two solutions have to be matched together, such 
that the curvature of JBJ for large arguments agrees with 
that of UlUfl for small arguments. At large capillary num- 
ber this will no longer be possible, since k££ becomes 
larger than fi9/9 e . At the critical capillary number Ca cr , 
© thus has to be equated with [J.9/9 e , giving 

SH 3 exp[-l/(35 cr )] = 2 1 /\(Ai(s max )) 2 f i6/9 e (11) 

for the critical rescaled capillary number S cr = 3Ca cr /9g. 
This equation completely determines the stability bound- 
ary found numerically in |l8j , as function of all parame- 
ters. It also implies that for a given geometry, instability 
occurs if the parameter S cr = ZCaj9\ exceeds a criti- 
cal value (up to very small logarithmic corrections), in 
agreement with experiment [l3j. 

For our procedure to be consistent, though, we need 
to make sure that the two solutions (JHJ and IjlOp have 
sufficient overlap to be matched. The inner and outer 
solutions at Ca cr are, in summary, 



H, 



Icr 

inner 
Icr 

outer 



= 5 1/3 /(^ 2 ) 

= (0/0 e )(l - e-*), 



where / is a universal function given by © with a = a cr . 
Thus ii P 2 ^> fj, the large-£ limit of H^~ ner overlaps with 
the small-£ limit of H'J£ ter , which translates into 9/9 e ^> 

s 1 / 3 . 

In Fig. [3we show the result of a numerical integration 
of for 9/9 cr = 10 at the critical capillary number. The 
interface slope follows the static solution perfectly almost 
up to the contact line, where it has to turn over, while 
the static solution extrapolates to zero. Coming from 
the interior, (JSJ describes the solution equally well up 
to the turning point. Note that there are no adjustable 
parameters in Fig. we simply took the inner and outer 
solutions in the critical case. In fact, even for 9/9 cr = 
1, when there is not yet much overlap between the two 
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FIG. 2: A comparison of the full solution at the critical 
capillary number with the inner and outer solutions. We plot 
the rescaled slope, so H'(0) = 1 for the full solution, and 
H'(0) = for the outer solution, consistent with the condition 
by Derjaguin and Levi. The parameters are fi = 10 -4 and 



= 10, and thus S cr = 0.0572 from JTTJk 



solutions, equation Ijlljl already works extremely well in 
describing the loss of the stationary solution, as shown 
in Fig. [3J Again, no parameter was adjusted to achieve 
this comparison. 

Beyond the wet-dry transition studied here, in two 
more states have been described, characterized by thin 
and thick wetting films, respectively. The stability of 
these solutions, and possible hysteresis phenomena re- 
main to be investigated. It is important to note that 
our approach is not limited to the moving plate geome- 
try studied here, nor is it restricted to a specific contact 
line model, since it is based entirely on hydrodynamic 
arguments away from the contact line. For example, if 
van-der-Waals forces are dominant near the contact line 
|2fij , the parameter fj, in (fTT|) simply needs to be replaced 
by Hvdw = \AV (67ry)/(20g), where A is the Hamakcr 
constant. 

To generalize to a different geometry, one has to replace 
(|10|l by the appropriate static solution for the problem 
at hand. This is done almost trivially for the case of a 
vertical wall [2l| , and easily extended to the flow in 
a narrow capillary, to be able to compare directly to ex- 
periments 01 . Hocking 0] found that the present slip 
theory correctly predicts Ca cr /9 3 to be a constant, but 
overestimates this constant by a factor of two, if reason- 
able values for the slip length A are assumed. Reasons 
for this discrepancy could be both the considerable con- 
tact angle hysteresis in the experiment |l3| . and failure 
of the simple contact line model used in our calculation. 
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FIG. 3: A comparison of S cr as determined by numerical 
integration of equation Q for 9 / 9 cr = 1 , and theory, as sum- 
marized by equation li lt . 



A new set of experiments, using the plate geometry, is 
being planned to clear up these questions |22j. 

Another important generalization is to higher dimen- 
sional problems, in which the contact line does not re- 
main straight. An instability toward inclined contact 
lines was observed in |12(, as well as in recent experi- 
ments with drops running down an inclined plane |23l |. 
To explain the characteristic inclination angle of such a 
contact line, one needs to identify a characteristic speed 
of de- wetting |12i 12 1| , which our present approach effec- 
tively provides. 

Special thanks are due to Petr Braun for numerous 
insightful discussions. 
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